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We investigate the gauge boson propagator in three dimensional compact Abelian gauge model 
in the Landau gauge at finite temperature. The presence of the monopole plasma in the confine¬ 
ment phase leads to appearance of an anomalous dimension in the momentum dependence of the 
propagator. The anomalous dimension as well as an appropriate ratio of photon wave function 
renormalization constants with and without monopoles are observed to be order parameters for 
the deconfinement phase transition. We discuss the relation between our results and the confining 
properties of the gluon propagator in non-Abelian gauge theories. 
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Three-dimensional compact electrodynamics (CQED 3 ) 
shares two outstanding features of QCD, confinement |f]] 
and chiral symmetry breaking |Q. With some care, 
it might be helpful for the understanding of certain 
non-perturbative aspects of QCD to study them within 
CQED 3 . The non-perturbative properties of CQED 3 
deserve interest by themselves because this model was 
shown to describe some features of Joseplison junc¬ 
tions ||] and high-T c superconductors 0|. 

Here, we want to elaborate on CQED 3 as a toy model of 
confinement. Indeed, this has been the first non-trivial 
case in which confinement of electrically charged particles 
was understood analytically [Q. Confinement is caused 
here by a plasma of monopoles which emerge due to the 
compactness of the gauge field. Other common features 
of the two theories are the existence of a mass gap and of 
a confinement-deconfinement phase transition at some 
non-zero temperature. According to universality argu¬ 
ments [|] the phase transition of CQED 3 is expected to 
be of Kosterlitz-Thouless type || . 

In QCD 4 , the deconfinement phase transition is widely 
believed to be caused by loss of monopole condensation 
(for a review see Ref. [ff ) within the effective dual su¬ 
perconductor approach Studying the dynamics of 
the monopole current inside gluodynamics, monopole de¬ 
condensation at the critical temperature is appearing as 
de-percolation, i.e. the decay of the infrared, percolat¬ 
ing monopole cluster into short monopole loops . This 
change of vacuum structure has a dimensionally reduced 
analog in the 3 D monopole-antimonopole pair binding 
which has been observed in cQED 3 [|o[ |n) . 

At present, the gluon propagator in QCD 4 is under 
intensive study. The analogies mentioned before encour¬ 
aged us to study the similarities between the gauge boson 


propagators in both theories. In order to fix the role of 
the monopole plasma in CQED 3 , not just for confinement 
of external charges but also for the non-perturbative 
modification of the gauge boson propagator, we consider 
it in the confinement and the deconfined phases. On the 
other hand, on the lattice at any temperature we are able 
to separate the monopole contribution to the propagator 
by means of eq. d) below. 

We have chosen the Landau gauge since it has been 
adopted in most of the investigations of the gauge boson 
propagators in QCD [|l2|, |l3| and QED [14, [lt|. In order 
to avoid the problem of Gribov copies 16], the alterna¬ 
tive Laplacian gauge has been used recently []l7j. The 
Coulomb gauge, augmented by a suitable global gauge 
in each time slice (minimal Coulomb gauge) has been 
advocated both analytically Q and numerically ©• 

The numerical lattice results for gluodynamics show 
that the propagator for all these gauges in momentum 
space is less singular than p ~ 2 in the immediate vicin¬ 
ity of p 2 = 0. Moreover, the results for the propagator 
at zero momentum are ranging from a finit e (lT j (Lapla¬ 
cian gauge) to a strictly vanishing |[g|, [h| [19| (Coulomb 
gauge) value. Recent investigations in the Landau gauge 
show that, beside the suppression at p —> 0 , the propaga¬ 
tor is enhanced at intermediate momenta which can be 
characterized by an anomalous dimension p~2| (see the 
last reference in Jl^] for a comparison of different model 
functions). 

In the present letter we demonstrate that the momen¬ 
tum behaviour of the photon propagator in QED 3 is also 
described by a Debye mass and by an anomalous dimen¬ 
sion which both vanish at the deconfinement transition. 
This mechanism can be clearly attributed to magnetic 
monopoles. The plasma contribution is relatively easy 
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to exhibit by explicit calculation and can be eliminated 
by monopole subtraction on the level of the gauge fields. 
The results of a study of the propagator in SU(2) glu- 
odynamics have been interpreted [Iij| ] in a similar spirit, 
where P-vortices appearing in the maximal center gauge 
were shown to be essential for the enhancement of the 
Landau gauge propagator at intermediate momenta. 

For our lattice study we have adopted the Wilson ac¬ 
tion, S[0] = f3J2p (1 — cos d p ), where 9 p is the U( 1) field 
strength tensor represented by the plaquette curl of the 
compact link field 9i, and /3 is the lattice coupling con¬ 
stant related to the lattice spacing a and the continuum 
coupling constant g 3 of the 3 D theory, (3 = l/(ag 3 ). We 
focus here on the difference between confined and decon- 
fined phase. All results presented have been obtained 
on lattices of size 32 2 x 8. The finite temperature phase 
transition is known to take place 00 at /3 C « 2.35. 

The Landau gauge fixing is defined by maximizing the 
functional ]F) ; cos0p over gauge transformations G. 
For details of the Monte Carlo algorithm we refer to 
[fiTf . A more complete presentation of our studies, in¬ 
cluding also a thorough analysis of the propagator in 
the zero temperature case is in preparation [plf . Details 
on the implementation of Landau gauge fixing, includ¬ 
ing the elimination of zero momentum modes and the 
careful control of double Dirac strings can be found in 
Ref. H f§. 

We study the gauge boson propagator, {9 p {x)9 l/ { 0)}, 
in the momentum space. The propagator is a function 
of the lattice momentum, p p = 2 sin(7rfc A1 /L /J ), where 
kp = 0,..., Lp /2 is an integer. We discuss here the finite 
temperature case and focus on the temporal component 
of the propagator, 

£>33(P 2 ,0) = —(0 3 (p,O)0 3 (-p,O)) (1) 

L/ x L J y± J z 

as function of the spatial momentum, p 2 = pfi ■ We 

remind that at finite temperature the confining proper¬ 
ties of static electrically charged particles are encoded in 
the temporal component of the gauge boson field, 9 3 . 

In order to pin down the effect of monopoles we have 
divided the gauge field 9i into a regular (photon) and 
a singular (monopole) part which can be done following 
Ref. [^2| . In the notation of lattice forms this is written: 

0 = 0P hot + d mon , 9 mon = 2ttA ~ l 5p{j) , (2) 

where A -1 is the inverse lattice Laplacian and the 0- 
form *j £ Z is nonvanishing on the sites of the dual lat¬ 
tice occupied by the monopoles. The 1-form *p[j} corre¬ 
sponds to the Dirac strings (living on the links of the dual 
lattice) which connect monopoles with anti-monopoles, 
S*p[j] = *j. For a Monte Carlo configuration, we have 
fixed the gauge, then located the Dirac strings, p[j] / 0, 
and constructed the monopole part 9 mon of the gauge 
field according to the last equation in (|^). The photon 


field is just the complement to the monopole part accord¬ 
ing to the first equation of (|j). 

The photon and monopole parts of the gauge field con¬ 
tribute to the propagator, D = D phot + D mon + D mlx , 
where D mlx represents the mixed contribution from reg¬ 
ular and singular fields. We show the propagator for 
p = (p, 0 ) together with the separate contributions, mul¬ 
tiplied by p 2 and averaged over the same p 2 values, in 
Figure [I] for coupling constant [3 = 1.8. 



FIG. 1: Different contributions to the full D 33 propagator 
(multiplied by p 2 ) vs spatial lattice momentum squared and 
fits as described in the text for /3 = 1.8 on a 32 2 x 8 lattice. 

The regular part of the propagator has perfectly the 
free field form 


7-vphot _ 

^33 — 


l ^Phot 


(3 p 2 


(3) 


at all available /?. The perturbative propagator defined in 
terms of 9i is obviously proportional to g 3l which is taken 
into account by the factor 1/(3 in eq. ([§[). The fits of the 
photon part of the propagator by the above expression 
give the parameter Z pbot as a function of lattice coupling 
(dash-dotted line in Figure [j]for (3 = 1.8). 

The singular contribution to the gauge boson propa¬ 
gator shows a maximum in p 2 D^ on at some momentum 
(Figure @), moving with increasing [3 nearer to |p| a = 0. 
The mixed component gives a negative contribution to 
p 2 D^ lx , growing with decreasing momentum. The cen¬ 
tral point of our paper is that all these contributions to¬ 
gether do not sum up to a simple massive Yukawa propa¬ 
gator. To quantify the difference between a Yukawa-type 
and the actual behavior we use the the following four- 
parameter model function for D 33 (p 2 ,0), 


d 33 (p 2 ,0) 


Z m 2a 

/3p2(l+“) +TO 2 (l+“) 


+ C, 


(4) 
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where Z , a, to and C are the fitting parameters. This 
model is similar to some of Refs. 12, [23| where the prop¬ 


agator in gluodynamics has been studied. 

The first part of the function (||) implies that the pho¬ 
ton acquires a Debye mass m (due to screening Q|) to¬ 
gether with the anomalous dimension a. The (squared) 
photon wave function renormalization constant Z de¬ 
scribes the renormalization of the photon wave function 
due to quantum corrections. The second part of ([ij) rep¬ 
resents a (5-function-like interaction in coordinate space. 

Before fitting we average the propagator over all lat¬ 
tice momenta at same p 2 to improve rotational invari¬ 
ance. Thus the errors entering the fits include both the 
variance among the averages for individual momenta and 
the individual errors. The fits were performed using stan¬ 
dard Mathematica packages combined with a search for 
the global minimum in \ 2 /d.o.f. To check the stability 
of the fits, we studied several possibilities of averaging 
and thinning out the data sets, a procedure which will 
be discussed elsewhere f2lj . 

The model function (|j) works perfectly for all p 2 and 
couplings (3. For /3 > 2.37 the best fit for mass parameter 
to and anomalous dimension a are both consistent with 
zero. Therefore we set to = 0 and a = 0 for these values 
of (3 to improve the quality of the fit of Z and C. 

It turns out that the inclusion of a constant term, (7, 
in the model function i) is crucial for obtaining good 
fits in the confinement phase, despite the fact that it is 
very small (as function of [3 the parameter C decreases 
from (7(1.0) = 0.18(4) to (7(2.2) = 0.009(2), it rapidly 
vanishes in the deconfined phase). Similarly to to and 
a parameters we set C to zero for (3 > 2.45, where C 
becomes smaller than 10 -4 . 

An example of the best fit of the full propagator for 
f3 = 1.8 is shown in Figure [l] by the solid line (with 
C = 0.033(5)). The parameter Z distinguishes clearly 
between the two phases (Figure [|). It coincides with 



FIG. 2: Coefficients Z of fit (|jj) for full propagator and Z phot 
for photon contribution (|j^) vs (3. 

the photon part Z phot (defined without monopoles) in 
the deconfined phase while it is much larger in the con¬ 
fined phase. This indicates that the photon wave function 


gets strongly renormalized by the monopole plasma. In 
contrast, the factor Z phot smoothly changes crossing the 
deconfinement transition at (3 C ~ 2.35. 

The anomalous dimension a also distinguishes the two 
phases (Figurep): it is equal to zero in the deconfinement 



FIG. 3: Anomalous dimension a vs (3 and its best fit near (3 C 
using function (jfij). 

phase (perturbative behaviour) while in the confinement 
phase the monopole plasma causes the anomalous dimen¬ 
sion growing to a ss 0.25 ... 0.3. 

To characterize the properties of Z and a approaching 
the phase transition we fit the excess of the ratio of Z’s 
over unity, 

= Z^W) ~ 1 ■ (5) 

and the anomalous dimension a in the following form: 

fi(P) = hi(P® -/?)*, (■i = a,Z ). ( 6 ) 

where i = Z,a. The f3i a ' Z ^ are the pseudo-critical cou¬ 
plings which might differ on finite lattices. 

The best fits f a and fz are shown in Figures || and 
[I], respectively. The solid lines in both plots extend over 



FIG. 4: Same as in Figure for ratio Rz, eq. (|E]). 

the fitting region. The corresponding parameters are pre¬ 
sented in Table |. The pseudo-critical couplings (3^ 
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i 

hi 

0 

7i 

a 

0.250(9) 

2.363(3) 

0.50(2) 

z 

2.63(7) 

2.368(5) 

0.48(3) 


TABLE I: Best parameters for the fits (^). 


(Z) 

and (3c ; are in agreement with previous numerical stud¬ 
ies 0 0 giving /3 C = 2.346(2). Note that the critical 
exponents 7 i are close to 1 / 2 , both for the anomalous 
dimension a and for Rz expressing the ratio of photon 
field renormalization constants. 

Finally, the /3-dependence of the mass parameter, m, 
is presented in Figure |j. As expected, the mass scale 
generated is non-vanishing in the confinement phase due 
to presence of the monopole plasma Q. It vanishes 
at the deconfinement transition point when the very 
dilute remaining monopoles and anti-monopoles form 
dipoles ]TT]. 



FIG. 5: The mass m vs /?■ 

Summarizing, we have shown that the presence of the 
monopole plasma leads to the appearance of a non¬ 
vanishing anomalous dimension a > 0 in the boson prop¬ 
agator of CQED 3 in the confinement phase. We would 
hope that our observation stimulates an analytical ex¬ 
planation. 

At this stage of studying cQED 3 as a model of con¬ 
finement we conjecture that in the case of QCD the 
Abelian monopoles defined within the Abelian projec¬ 
tion may be responsible for the anomalous dimension of 
the gluon propagator observed in Refs. jl^, ^3). If true, 
a monopole subtraction procedure analogous to that em¬ 
ployed here would be able to demonstrate this. We found 
that the anomalous dimension a and the ratio of the pho¬ 
ton wave function renormalization constants with and 
without monopoles, Rz (^j), represent alternative, also 
non-local order parameters characterizing the confine¬ 
ment phase. 
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